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ABSTRACT 


The  electromechanical  behavior  of  a  Single 
Stage  Induction  Mass  Driver  (SSIMD)  has  been  modeled 
and  analyzed  with  the  aid  of  a  computer .  The 
analysis  focused  on  the  effects  of  circuit 
parameters  on  projectile  velocity  and  system 
efficiency.  A  description  of  the  model  along  with 
general  rules  for  the  design  of  efficient  SSMID 
systems  are  presented. 
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INTRODUCTION 


The  promise  of  achieving  high  velocities  and 
high  efficiencies  has  made  the  SSIMD  an  attractive 
concept  for  many  applications.  Since  the  mid¬ 
sixties  the  SSIMD  has  been  adapted  to  such 
applications  as:  magnetic  flux  compressors  {l}, 
hypervelocity  mass  accelerators  {2},  induction 
reaction  engines  {3}, {4},  and  electromagnetic  metal 
formers  {5} . 
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Figure  1 :  Single  Stage  Induction  Mass  Driver 


The  generic  SSIMD  (figure  1)  basically  consists 
of  an  energy  storage  (normally  a  capacitor  bank),  a 
switch,  a  drive  coil,  and  a  driven  ring.  The  coil 
and  ring  are  concentric  and  typically  of  equal 
diameter  ( the  condition  for  maximum  magnetic 
coupling) .  When  a  current  pulse  is  applied  to  the 
drive-coil,  a  current  is  magnetically  induced  in  the 
ring.  The  interaction  between  this  ring  current  and 
the  magnetic  field  (between  the  coil  and  the  ring) 
produces  a  Lorentz  Force  which  drives  the  ring. 

The  ring  and  coil  geometries,  and  their 
relative  orientation  vary  depending  on  the  desired 
drive-direction  (figure  2).  In  figure  2a  the  ring 
is  driven  axially,  away  from  the  drive  coil.  This 
is  the  geometry  normally  used  in  reaction  engines 
and  hypervelocity  accelerators.  In  figure  2b  the 
ring  is  driven  inward  radially,  making  this  coil¬ 
ring  arrangement  favorable  for  flux-compression  and 
metal  forming  applications. 

In  an  effort  to  explore  efficient  induction 
driving,  a  twin  capacitor  bank  system  has  been  built 
and  computer  simulated.  This  twin  bank  scheme 
allows  the  use  of  high  energy-density  capacitors 
(for  a  compact  system)  by  preventing  voltage 
reversal  across  the  banks.  This  paper  is  based 
mainly  on  experimental  and  simulation  work  with  the 
above  system. 


Figure  2:  The  two  basic  drive  system  geometries, 

a.  typical  geometry  of  axial  accelerators 

b.  typical  geometry  of  flux  compressors 
and  metal  formers 

CIRCUIT  MODEL 


Our  system  has  been  electrically  modeled  by  the 
circuit  shown  in  figure  3.  For  discussion  purposes, 
the  circuit  has  been  divided  into  two  main  sections: 
( 1 )  the  drive  system,  which  includes  the  drive  coil 
and  the  ring;  and  (2)  the  power  system,  which 
includes  the  capacitor  banks,  the  switch  and  the 
power  bus. 


Figure  3 :  Circuit  model  of  twin  capacitor-bank 
device 

The  drive  system  has  been  modeled  as  a 
transformer  with  a  shorted  secondary  and  a  mutual 
inductance  (M)  which  is  a  function  of  ring 
separation.  The  self  inductance  of  the  primary  (Li) 
is  the  drive  coil  Inductance.  This  value  is 
proportional  to  the  coil  mean  diameter  and  the 
number  of  turns  squared.  The  secondary  self 
inductance  (L2)  is  the  ring  inductance,  which, 
within  a  few  percent  can  be  approximated  by  {6} 

r  ,Rm.  1  (1) 

L  =y  R  -Un(xj-)+.9} 

2  o  m  w 

Where  R^  is  the  ring's  mean  radius  and  w  is  its 
width.  L2  is  a  variable  quantity  if  ring  deformation 
occurs  during  acceleration.  Due  to  the  typically 
high  turns  ratio  between  the  primary  and  the 
secondary,  the  ring  sees  the  highest  currents  in  the 
system,  and  therefore  experiences  significant 
heating.  For  this  reason,  R2,  the  ring  resistance, 
is  modeled  as  a  temperature  dependent  element 
defined  by  the  following  equation: 

R2=R2o{l+VT-To)}  (2) 
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Where  Kt  is  the  temperature  coefficient  of  the  ring- 
material  (for  copper  and  aluminum,  at  20  °C, 

Kt=-0039  {7}),  TQ  is  the  initial  ring  temperature,  T 


is  the  operating  temperature, 
ring  resistance  -defined  by 

2irR 


and  R20  is  the  initial 
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2o 


CwT, 
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where  C  is  the  ring-conductivity,  and  Tj,  is  the  ring 
thickness.  The  separation  dependent  mutual 
inductance,  M,  is  determined  using  the  following 
semiempirical  procedure:  First,  the  drive  system 
circuit  is  transformed,  ignoring  R2,  to  its  T-model 
equivalent  (figure  4). 


Figure  4 :  T-model  of  drive  system  used  to  obtain 
M(z) 


From  this  circuit,  the  effective  drive  system 
inductance,  Leff,  (measured  at  the  primary 
terminals)  is 


wv<$ 


(4) 


there  is  significant  coil-ring  magnetic  interaction, 
and  it  is  defined  as  the  separation  when  Lef  .9L-| . 
This  coupling  length  has  experimentally  been  found 
to  roughly  equal  .17D,  where  D  is  the  mean  drive- 
coil  diameter.  This  result  agrees  with  streak 
camera  observations  by  Bondaletov  {2}  of  the  maximum 
separation  up  to  which  the  ring  is  accelerated. 

In  the  power-system,  R-j  is  the  lumped 
representation  of  the  primary  loop  resistance,  and 
it  is  obtained  from  the  exponential  decay  of  the 
oscillating  primary  current  (1-j)  without  the  ring  in 
place.  Similarly,  the  stray  inductance,  Ls,  is 
obtained  from  the  oscillating  frequency  of  1^  when 
the  coil  is  replaced  by  a  shorting  strap.  Finally, 
C,  and  C2  represent  two  capacitor  banks  of  equal 
capacitance;  where  initially  C-|  is  charged  and  C2  is 
uncharged. 


SIMULATION 


In  order  to  predict  the  behavior  of  our  SSIMD, 
a  computer  simulation  has  been  developed  using  the 
circuit  model  previously  described.  The  time 
behavior  of  the  following  variables  were  of 
interest:  the  primary  current  (I-]),  the  secondary 
(ring)  current  (I2),  the  capacitor  banks  voltages 
(Vi  and  V2).  the  ring's  normal  velocity  (V ) ,  the 
normal  force  exerted  on  the  ring  (F),  the  normal 
ring  position  (z),  and  the  heat  energy  dissipated  in 
the  ring  (E2).  The  approach  taken  was  to  first 
write  differential  loop  equations  for  the  primary 
and  the  secondary  of  figure  3,  and  then  manipulate 
them  into  a  form  that  could  be  solved  numerically. 
For  the  primary 


(5) 


Llil+IlR1^2/C2"VCl+MVVI=0 


By  measuring  Leff  at  increments  of  ring-coil  normal 
separation,  z,  curves  for  Leff  and  M,  similar  to  the 
ones  in  figure  5,  are  obtained. 


Figure  5:  Experimental  curves  for  drive  system 
effective  inductance  and  mutual 
inductance  (coil  mean  diameter 
=20.5  cm,  6  turns) 

Finally,  an  exponential  curve  fit  to  the  discrete 
M(z)  and  beff(z)  yields  the  following  expressions 


L  f=L(l-Aa2e-2z/2t] 
eff  1  L2 


-zAo 


L.  L_ 
2o  1 


K  e 
o 


(6) 


(7) 


where  K0  is  the  initial  (at  z=0)  coupling 
coefficient  (as  used  in  transformer  terminology)  and 
is  defined  by 


K  =/l-C 
o  f 


(8) 


and  the  secondary 

(i2+R2)I2+L2i2+Mil+IlA=0 


(10) 


Where  Qi  and  Q2  are  respectively  the  charge  stored 
in  Ci  and  C2.  The  numerical  technique  used  was  a 
fourth  order  Runge-Kuta  Method,  which  requires  that 
the  derivative  of  the  variable  to  be  solved  be  equal 
to  an  expression  which  does  not  include  that 
derivative.  By  respectively  multiplying  equations 
(9)  and  (10)  by  L2  and  M,  then  subtracting  the 
resultant  equations  an  expression  for  dl-|/dt  is 
obtained 


WVWV  (Mft-L2Ri]  1 

VrM 

+(ml2+r2m-l2m)i2> 


1 


(11) 


Similarly,  by  respectively  multiplying  equations  (9) 
and  (10)  by  M  and  L^,  and  then  subtracting  the 
resultant  equations  dl2/dt  is  obtained 


i2= - ^-2-{mq2/c2-mq1/c1+(mr1-ml1)  i1 


L2Li-M 


(12) 


+(mk-l1L2-r2l1)i2} 

The  rest  of  the  necessary  equations  are : 


(  13) 


where 


V1=WC1 


where  Cf,  the  compensation  factor,  is  the  factor  by 
which  L-j  is  reduced  when  the  ring  is  at  z=0.  Z^, 

the  coupling  length,  is  the  separation  up  to  which 
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(14) 


where 


70 


WC2  (16) 

The  force  exerted  on  the  ring  is  defined  by  {-8 } 

F=I1I2dK/dz  (17) 


or 


V  =  F/m- (I  I2/m) (dM/dz) 

where  m  is  the  mass  of  the  ring. 

z=V 


(18) 


(19) 


and 

fe2=I22R2  <20> 

Figure  6  shows  the  various  simulation  output  curves 
for  a  sample  run  along  with  comparable  measured 
variables. 


60 


Figure  7 :  Computer  generated  optimum  efficiency 

curves  as  a  function  of  Li,  Ceff,  and  Ls. 
Rl-15.6mft,  R2-1.1m0,  L2=.32iiH,  ^=.89, 
Z£=. 035m,  and  m=.1  Kg 


Figure  6s  Computer  generated  curves  compare  to 
measurements  for  the  following  case: 
Li"14)iH,  Ci~C2=828)IF,  Ri~34mfi,  Lg-lpH, 
2^».035m,  L2o=.32JiH,R2o».65mfi,  m=.115  Kg, 
and  Kq“ .89 

In  general,  good  agreement  exists  between  measured 
and  calculated  values.  The  most  signficant  error 
exists  in  the  velocity  prediction  (typically  by  a 
factor  of  ~ . 7 ) .  This  error  is  probably  due  to 
aerodynamic  effects  which  have  not  been  modeled. 


ANALYSIS 


In  an  effort  to  emperically  develop  rules  for 
the  design  of  energy  efficient  SSIMD,  the  previously 
described  simulation  was  used  to  generate  the  family 
of  curves  shown  in  figure  7 . 


The  solid  and  dashed  curves  respectively  represent 
the  presence  and  absence  of  stray  inductance  in  the 
system,  and  have  been  generated  for  three  values  of 
Ci  (C2«Ci>. 

Each  point  in  these  curves  represents  an 
optimum  efficiency  for  a  given  set  of  circuit 
elements,  and  was  found  by  stepping  through  the 
charge  voltage  until  an  efficiency  maximum  was 
reached.  The  purpose  of  changing  the  voltage  was  to 
vary  the  ring  velocity.  These  points  agree  with  the 
following  criterion  { 4) : 

<V>— z^/<T/2>  (21) 

where  <\»  is  the  average  ring  velocity  and  <T/2>  is 
the  average  first  half-period  of  the  primary  current 
-  defined  by 


<T/2>=7r/<Leff>Ceff 


(22) 


where 


3ff>: 


Jo 


L  „dz+L 
ef  f  s 


(23) 


and  Ceff  is  the  effective  system  capacitance.  For 
example,  in  the  twin  capacitor  bank  system  Ceff 
equals  the  series  combination  of  Ci  and  C2.  The 
above  criterion  simply  states  that  in  order  to  drive 
a  ring  most  efficiently,  the  ring  must  remain  in  the 
coupling  region  (between  z=0  and  z=Z^ )  for  exactly 


<  T/2  >  . 


By  comparing  the  solid  and  dashed  curves  the 
adverse  effects  of  stray  inductance  on  efficiency  is 
clearly  seen.  When  Ls  is  of  the  order  of  Li  a 
signficant  energy  division  occurs  between  the  two: 
this  division  reduces  the  amount  of  energy  available 
at  the  drive  coil,  and  consequently  reduces  drive 
efficiency.  Aside  from  the  effects  of  Ls  on 
efficiency,  when  Li  is  made  too  small,  the 
efficiency  is  also  reduced,  as  shown  by  the  dashed 
curves  of  figure  7.  This  effect  can  be  attributed 
to  the  effective  coupling  parameter  dM/dz |  eff,  which 
is  present  in  the  force  equation,  as  shown  below 


F= 1 1 12dM/dz=I 12M/L2dM/d2=I 12dM/dz  j 


eff 


(24) 


This  effective  coupling  parameter  is  used  as  a 
figure  of  merit  when  comparing  difference  induction 
accelerators  { 4  }.  The  dependence  of  dM/dz | ef f  on  Li 
can  be  seen  when  M  and  dM/dz  are  spacially  averaged 
from,  z=0  to  and  substituted  into  eqn.  (24) 


>=- 

eff 


.  4K  2L 
o  1 


(25) 


although  it  is  undesirable  for  Li  to  be  small,  it  is 
also  undesirable  for  it  to  be  too  large,  since  Ii, 
(the  other  term  in  eqn.  (24)),  will  be  reduced. 

These  curves  also  show  the  beneficial  effects, 
in  terms  of  efficiency,  of  reducing  the  effective 
system  capacitance.  Since  lower  capacitance 
translates  to  higher  velocity  (from  the  optimum 
efficiency  criteria),  the  trends  in  figure  7  imply 
that  induction  drivers  have  a  natural  tendency  to 
efficiently  drive  rings  to  high  velocities. 
Unfortunately,  the  price  paid  for  reducing  the 
capacitance  and  still  operate  at  optimum  efficiency 
is  an  increase  in  charge  voltage,  which  ultimately 
becomes  one  of  the  practical  limits  on  device 
performance.  This  voltage  limit  is  approximately 
lOOKv  {4  }.  An  estimate  of  the  optimum  charge 
voltage,  which  agrees  with  the  simulation  results, 
can  be  derived  from  the  force  equation  (eqn.  (17)) 
by  averaging  all  of  the  pertinent  variables:  if 
L2/R2  »<T/2>,  which  usually  is  the  case,  eqn.  (17) 
can  be  rewritten  as 
2 

F=I^  (M/L2)dM/dz=ma  (26) 

Averaging  this  expression  and  replacing  <a>  ,  the 
average  ring  acceleration,  with  <v>/<T/2> 


(2m<V>)/<T/2>=I  /4<M/L_ (dM/dz) > 

P  2 


(27) 


where  Ii  was  replaced  by  its  RMS  value  from  t-0  to 
t=<T/2>,  Ip/1/4.  Ip  is  the  peak  primary  current, 
defined  by 


I  =' 
P 


1 


eff 


<L 

eff 


(28) 


Finally,  combining  eqns.  (25)  and  (26)  yields  an 
expression  for  the  optimum  charge  voltage. 


V  ■ 
o 
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4m<\» 


/<L  > 

eff 


eff  L2 


dM 

dz 


(29) 


The  curves  in  figure  7  are  independent  of  ring 
mass,  but  the  charge  voltage  required  to  operate  at 
optimum  efficiency  is  proportional  to  the  square 
root  of  the  mass  -as  indicated  by  Eqn.  (29). 

DESIGN  PROCEDURE 


The  purpose  of  this  procedure  is  to  determine 
the  values  of  the  circuit  elements  and  the  charge 
voltage,  based  on  a  design  velocity  and  drive  system 


dimensions,  wnich  will  yield  an  efficient  SSIMD. 

The  objective  is  to  choose  the  circuit  parameters 
which  will  cause  the  device  to  operate  near  the 
optimum  efficiency  peaks  of  figure  7.  The  procedure 
assumes  a  capacitive  discharge,  series  RLC-circuit ; 
and  a  drive-coil  and  ring  of  roughly  equal  diameter. 
The  first  step  is  to  minimize  the  values  of  Ls  and 
Rl .  Lg  is  determined  by  the  geometry  of  the  power 
system;  some  typical  low  values  for  parallel  plate 
and  coaxial  geometries  are  respectively  250nH{3}  and 
120nH  4  (at  voltages  near  30Kv).  Ri  is  the 
conductor  resistance  of  the  primary  circuit,  and  is 
typically  of  the  order  of  10E-3  ohms.  The 
application  normally  dictates  the  dimensions  of  the 
drive-coil  and  the  ring.  These  dimensions  in  turn 
determine  Lj  and  the  lower  limit  of  Li .  By 
increasing  the  number  of  turns,  Li  is  increased 
until  Li  > 10LS.  The  upper  limit  on  Li  is  not  easily 
defined,  and  it  may  be  necessary  to  use  a  computer 
model  to  determine  its  optimum  value.  R2  is 
determined  by  the  ring's  dimensions  and  material 
properties.  This  parameter  has  a  strong  impact  on 
efficiency  through  the  l22R2  loses,  and  therefore 
should  be  minimized.  Ceff  is  chosen  according  to 
the  optimum  efficiency  criteria  (eqn.  (21)),  where 
the  average  velocity  can  be  approximated  by  one  half 
the  design  velocity.  Finally,  a  rough  estimate  of 
the  required  voltage  is  obtained  using  eqn.  (29), 
where  the  spacial  averages  for  M,  Lef f ,  and  dM/dz 
are  found  from  their  previously  described  equations. 
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